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ABSTRACT 

A  conditional  moment  closure  model  is  proposed  for  reacting  particles  in 
turbulent  nonpremixed  combustion.  The  new  model  for  particles  differs  signif¬ 
icantly  from  the  traditional  uniform  diffusivity  gas-phase  conditional  moment 
closure  model.  The  new  features  of  the  model  and  model  effectiveness  are 
examined  against  direct  numerical  simulation  data  for  soot-like  and  droplet- 
like  particles  in  turbulent  nonpremixed  combustion.  The  influence  of  differing 
particle  sizes,  types,  and  reactivity  on  the  effectiveness  of  the  model  closure  is 
examined  in  detail. 
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Simulation  and  modelling  of  reacting  particles  relevant  to 
gas  turbine  engine  combustion 


EXECUTIVE  SUMMARY 

The  purpose  of  DST  Task  98/091  Modelling  Airbreathing  Combustion  has  been  to 
develop  a  capability  to  numerically  predict  important  characteristics  in  generic  gas  tur¬ 
bine  combustion  systems  using  a  fundamentally  sound,  and  therefore  generally  applicable, 
method. 

Liquid  fuel  droplets  and  solid  carbonaceous  soot  particles  are  two  key  particle  types 
which  frequently  occur  in  gas  turbine  engine  combustors  during  operation.  Fuel  droplets 
are  evaporated  by  the  transfer  of  heat  to  the  liquid  from  hot  combustion  gases  and  heated 
surfaces.  The  distribution  of  combustible  gaseous  fuel  is  strongly  dependent  on  the  dis¬ 
tribution  of  fuel  droplets  and  their  rate  of  evaporation.  Poor  fuel  evaporation  can  lead 
to  poor  combustion  stability  and  efficiency,  and  have  a  strong  detrimental  influence  upon 
engine  operation. 

The  formation  of  soot  and  its  impact  upon  radiation  heat  transfer  from  combustion 
systems  and  their  hot  exhausts  is  of  critical  significance  to  the  maintenance  and  tactical 
effectiveness  of  modern  military  aircraft.  Radiation  from  soot  is  the  key  heat  transfer 
mechanism  from  combusting  gases  to  surrounding  solid  surfaces  in  the  hot  section  of 
gas  turbine  engines.  The  level  and  distribution  of  this  heat  transfer  plays  a  determining 
role  in  the  life  of  hot  section  components,  most  notably  the  combustor  liner.  Further, 
the  presence  of  hot  soot  particles  in  the  exhaust  stream  from  an  aircraft  powerplant  is 
thought  to  influence  the  nature  of  the  infra-red  signature  of  the  aircraft,  perhaps  adding 
a  substantial  grey  body  emission  source  to  the  exhaust  plume  itself. 

To  date,  soot  modelling  has  been  the  primary  focus  of  task  98/091,  as  a  result  of 
its  importance  to  predicting  engine  operation  characteristics,  as  described  above.  In  a 
prior  report,  entitled  Comparison  of  soot  model  predictions  with  experimental  data  for  a 
turbulent  sooting  propane  jet  flame  (DSTO-TR-0676),  the  first  efforts  at  modelling  sooting 
combustion  using  a  Conditional  Moment  Closure  (CMC)  method  were  described.  In  that 
report,  experimental  data  for  turbulent  sooting  propane  jet  flames  were  compared  with 
CMC  model  predictions.  While  the  agreement  between  the  model  and  experimental  data 
was  quite  good  it  was  clear  that  a  better  understanding  of  soot  particle  modelling,  and  a 
better  means  of  treating  radiation  heat  transfer,  was  required. 

This  report  discusses  the  derivation  and  validation  of  a  modified  CMC  model  which  is 
designed  to  account  for  the  reactions  of  particles,  such  as  soot  and  droplets,  in  turbulent 
nonpremixed  combustion.  Initial  validation  of  the  CMC-particle  model  was  achieved  by 
comparing  model  predictions  with  Direct  Numerical  Simulation  (DNS)  data  for  selected 
flow  and  mixing  conditions  which  represent  exact  numerical  solutions  for  idealised  cases 
of  turbulent  combustion.  It  was  found  that  the  model  has  significant  utility  in  predicting 
the  behaviour  of  reacting  particles  in  turbulent  combustion,  and  that  further  research  and 
development  in  this  area  is  warranted.  It  was  found  that  the  modelling  of  large  inertial 
particles  is  made  somewhat  problematic,  in  the  framework  of  the  CMC  model,  by  the 
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failure  of  these  particles  to  closely  follow  the  fluid  flow.  However,  flow-following  particles 
such  as  soot  and  smaller  droplets  were  modelled  very  successfully. 
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1  Introduction 

Condensed-phase  particles  are  frequently  present  in  turbulent  combustion  systems  and 
can  have  a  profound  influence  on  the  thermochemical  nature  of  their  surroundings.  Fuel 
droplets  and  soot  particles  are  two  examples  of  important  condensed  phase  species  in 
combustion.  The  evaporation  of  the  former  largely  determines  the  distribution  of  the 
combustible  gaseous  fuel/air  mixture,  while  the  presence  of  the  latter  impacts  strongly 
on  the  degree  of  radiant  heat  transfer  from  the  system.  It  is  thus  desirable  to  be  able  to 
predict  the  mean  behaviour  of  these  particles,  in  a  turbulent  combusting  environment,  in 
response  to  their  local  thermochemical  conditions. 

A  number  of  difficulties  surround  the  modelling  of  the  mean  rate  of  particle  reactions 
in  turbulent  combustion.  Perhaps  the  most  significant  difficulty  is  associated  with  deter¬ 
mining  the  mean  influence  of  carrier  fluid  properties  such  as  temperature  and  chemical 
species  concentrations  upon  the  particle  population.  Wherever  the  local  properties  of  the 
carrier  fluid  fluctuate  due  to  turbulence,  the  variations  can  couple  with  the  non-linear 
particle  reactions  to  preclude  first  order  closure  of  the  mean  rates  with  mean  properties. 
This  type  of  closure  problem  is  the  same  as  that  encountered  in  modelling  the  mean  rate 
of  purely  gas  phase  chemical  reactions  using  conventional  averaging  techniques. 

The  Conditional  Moment  Closure  (CMC)  method  (see  Klimenko[l],  Bilger[2])  for  mod¬ 
elling  turbulent  gas  phase  nonpremixed  combustion  makes  use  of  averages  which  are  con¬ 
ditional  upon  the  local  value  of  a  conserved  scalar  (mixture  fraction),  which  is  indicative 
of  the  state  of  mixing  between  fuel  and  air  masses.  Conditional  averaging  upon  mixture 
fraction  captures  much  of  the  turbulence-induced  fluctuations,  and  a  first  order  condition¬ 
ally  averaged  closure  is  often  possible.  The  success  of  the  CMC  model  in  predicting  gas 
phase  combustion  makes  it  of  some  interest  in  modelling  particle  reactions  in  turbulent 
flow. 

The  purpose  of  this  study  was  to  simulate  the  dynamics  of  reactive  particle  mass  and 
motion  in  a  turbulent  combusting  environment  in  an  attempt  to  model  the  observed  mean 
thermochemical  behaviour  of  the  reacting  particles  using  a  derivative  of  the  CMC  method. 
In  this  report,  the  simulation  and  modelling  of  pseudo-soot  and  pseudo-droplet  particles 
are  described.  The  former  type  of  particles  were  smaller  than  the  latter,  and  they  were 
subject  to  relatively  strong  processes  of  growth  and  consumption.  These  processes  were 
meant  to  represent  soot  surface  growth  and  oxidation.  The  pseudo-droplet  particles  were 
subject  to  an  evaporative  process  only.  In  both  the  simulations  and  modelling,  particle- 
particle  interactions  were  disallowed,  thus  rendering  the  results  of  this  study  valid  only  for 
low  particle  mass  loadings  of  the  gas  phase. 


2  Simulation  Conditions 

Direct  numerical  simulations  were  conducted  using  a  pseudo-spectral  solution  tech¬ 
nique  for  forced  isotropic  turbulent  flow  on  a  32-cubed  grid.  The  flow  was  incompressible 
in  nature,  but  a  passive  scalar  (mixture  fraction)  field  was  used,  in  conjunction  with  an 
equilibrium  temperature  profile  in  passive  scalar  space,  to  determine  a  false  temperature 
map  throughout  the  domain.  This  false  temperature  was  computed  for  the  purpose  of 
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determining  instantaneous  particle  reaction  rates  at  every  step.  The  passive  scalar  field 
was  forced  at  large  scales  through  the  interaction  of  turbulent  motions  with  an  imposed 
mean  scalar  gradient  in  the  x  direction,  in  the  same  manner  as  the  simulations  of  Over¬ 
holt  and  Pope[3].  Cubic  spline  tensor-products,  as  described  by  Yeung  and  Pope[4,  5], 
were  employed  to  determine  local  fluid  velocities  and  temperatures  at  all  particle  loca¬ 
tions.  Particle  properties  were  advanced  in  time  using  a  fourth-order  Adams-Bashforth 
timestepping  routine. 

The  following  Lagrangian  equation  for  particle  motion  was  solved, 

~^  =  ~  (™T2/3  («»(*)  -  v{)  ,  (1) 

where  u,  denotes  the  particle  velocity,  U{(x)  denotes  the  local  fluid  velocity  at  the  particle 
location,  m*  is  the  nondimensional  particle  mass,  r*  is  the  Kolmogorov  timescale,  and  a* 
is  the  nondimensional  characteristic  particle  rate  given  by, 


In  the  above,  Lfc  is  the  Kolmogorov  length  scale,  d  is  the  reference  particle  diameter,  pp  is 
the  particle  material  density,  and  pj  is  the  fluid  density. 

The  Stokes  drag  expression  above  (Eqn  1)  strictly  applies  only  in  the  limit  of  a  purely 
laminar  flow  around  the  particle.  This  assumption  is  valid  for  the  smaller  soot  particles 
but  is  not  as  well  justified  for  the  droplets.  Empirical  corrections  for  changes  in  particle 
drag  with  higher  slip  Reynolds  numbers  are  available  but  have  not  been  applied  in  the 
simulations  reported  here. 

For  droplets,  the  rate  of  change  of  the  mass  of  any  particle  was  given  by 


dm" 

dt 


(3) 


where  /?  is  a  generic  rate  coefficient,  nominally  equal  to  a*  for  inertial  particles,  and  T* 
is  a  nondimensional  temperature  equivalent  to  the  local  Spalding  transfer  number[ 6].  The 
value  of  the  temperature  varied  from  zero  in  non-reactive  regions  of  passive  scalar  space 
to  a  peak  value  of  3.4  at  stoichiometric  conditions.  The  peak  value  corresponds  to  the 
transfer  number  for  kerosene  droplets  evaporating  within  enveloping  flames  (see  Kuo[6j). 

The  rate  of  change  of  mass  for  any  given  soot  particle  was  given  by 


dm* 

dt 


(4) 


where  (5  is  a  rate  coefficient,  not  directly  related  to  a*,  and  fgr  and  fco  are  normalised 
functions  of  local  mixture  fraction  (and  thus  an  inferred  reactive  gas  phase  composition) 
which  mimic  surface  growth  and  consumption  respectively.  The  growth  and  consumption 
reactions  were  designed  to  embody  the  basic  features  of  the  soot  processes  they  represent, 
but  with  one  notable  difference.  Namely,  that  the  reaction  rate  profiles  were  defined  so  as  to 
be  symmetric  in  mixture  fraction  space  about  the  mean  mixture  fraction.  Particle  surface 
growth  was  strongest  at  the  mean  mixture  fraction,  and  particle  consumption  was  strongest 
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Case 

Type 

a* 

0 

Tsim/  ^eddy 

T sim / Tpop 

dl 

droplet 

0.014 

0.014 

193 

34.7 

dli 

droplet 

oo 

0.014 

96.1 

15.0 

si 

soot 

0.7 

2.8 

3.5 

0.46 

s2 

soot 

0.7 

0.7 

104 

23.1 

s2i 

soot 

00 

0.7 

104 

23.1 

Table  1:  Parameters  and  durations  for  the  reactive  particle  simulations. 

somewhat  to  the  lean  and  rich  sides  of  the  mean  mixture  fraction.  This  symmetry  of 
thermochemical  properties  about  the  mean  mixture  fraction  allowed  conditional  statistics 
on  either  side  of  the  mean  to  be  combined  to  increase  statistical  significance  in  each  half¬ 
plane  of  mixture  fraction  space. 

The  reference  particle  rates  (a*,  /3)  differed  between  simulations  to  reflect  different 
sizes  and  reactivity  of  the  particles,  whilst  flow  and  mixing  conditions  were  the  same  for 
all  cases.  The  details  of  the  simulation  cases  studied  are  given  in  Table  I.  Some  inertial 
particle  simulations  were  repeated  with  inertialess  particles  (denoted  by  i-suffixed  case 
designations  in  Table  I)  to  ascertain  the  influence  of  particle  slip  velocity  on  the  reactive 
particle  statistics. 

In  the  table,  the  ratios  of  the  duration  (rs,m)  of  each  case  to  both  the  large  eddy 
turnover  time  (reddy),  and  the  mean  time  for  total  particle  population  change  through 
injection  (rpop),  are  given.  Owing  to  explicit  method  time  step  limitations,  the  cases  with 
longer  characteristic  timescales,  such  as  the  droplet  and  less-reactive  soot  cases,  were  less 
expensive  to  compute.  These  cases  could  therefore  could  be  continued  for  a  larger  number 
of  eddy-turnover  times  for  a  given  amount  of  real  time. 

In  each  case,  8192  particles  were  maintained  in  the  domain,  at  a  constant  mean  number 
density,  by  continuous  injection.  Particle  simulations  commenced  only  after  the  forcing 
had  produced  a  statistically  stationary  flow'  field  over  a  period  of  many  hundreds  of  eddy 
turnover  times.  The  simulations  were  conducted  at  a  Taylor  Reynolds  number  of  ~  24, 
with  a  value  of  approximately  0.22  for  the  ratio  of  integral  length  scale  to  domain  width, 
and  2.4  for  the  product  of  Kolmogorov  lengthscale  and  the  maximum  wavenumber  in  the 
simulation. 

The  actual  size  of  the  soot  particles  used  in  the  inertial  simulations  of  this  study 
were  at  the  upper  limit  of  realistic  soot  particles.  In  a  sense,  these  simulation  cases 
embody  a  worst  case  scenario  for  modelling  soot.  In  practical  applications,  the  bulk  of  soot 
particles  are  much  smaller  and  follow  the  flow  field  to  much  more  closely.  These  smaller 
reacting  particles  were  better  simulated  using  an  inertia-less  particle  tracking  method  in 
conjunction  with  the  direct  numerical  simulations  of  turbulence. 


2.1  Particle  Injection 


The  imposed  mean  scalar  gradient  in  the  x  direction  of  the  simulation  domain  required 
that  special  consideration  to  be  given  to  particles  which  crossed  the  domain  boundaries 
in  this  direction.  In  the  other  directions,  all  properties  were  periodic  and  particles  were 
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simply  returned  to  the  domain  by  a  periodic  mapping.  The  same  mapping  was  not  applied 
in  the  x  direction  as  that  would  allow  particles  to  transfer  between  very  lean  and  very  rich 
mixture  conditions  in  crossing  the  domain  boundary. 

Instead,  any  particle  which  crossed  an  x  boundary  was  deemed  to  have  left  the  domain 
permanently,  and  was  replaced  via  a  periodic  mapping  by  a  fresh  particle  with  a  velocity 
equal  to  the  departing  particle,  and  mass  equal  to  the  injection  mass.  This  type  of  fresh 
particle  injection  allowed  a  statistically  stationary  particle  state  to  be  reached  as  a  bal¬ 
ance  between  the  processes  of  particle  reaction,  and  transport  to  and  from  the  injection 
boundaries  occurs. 

Injection  across  the  two  ^-bounding  planes  was  characterised  by  a  relatively  weak  in¬ 
fluence  on  conditional  mean  particle  statistics  at  very  rich  and  lean  mixture  fractions.  The 
source  terms,  for  the  conditional  mean  and  variance  equations  for  particle  mass  fraction, 
which  result  from  these  injection  schemes  are  described  in  Section  3.2. 

2.2  Mixing  and  Reaction  Statistics 

Statistical  data  based  on  particle  mass  fractions  were  measured  from  the  simulations 
for  comparison  with  model  predictions.  These  data  were  recorded  only  after  the  par¬ 
ticle  ensemble  had  been  determined  to  be  statistically  stationary,  by  monitoring  initial 
transients  in  the  data. 

As  each  particle,  on  average,  occupied  four  grid  cells  to  the  exclusion  of  other  particles, 
it  was  not  possible  to  determine  meaningful  spatial  gradients  in  particle  properties  within 
the  carrier  fluid.  It  was  possible,  however,  to  derive  conditional  statistical  information 
from  the  data  by  binning  particles  in  mixture  fraction  space  according  to  their  local  value 
of  mixture  fraction.  The  particle  density  in  each  bin  allowed  gradients  of  conditional 
quantities  in  mixture  fraction  space  to  be  determined  for  comparison  with  modelled  terms. 

Mixing  statistics  such  as  conditional  mean  scalar  dissipation  rate,  conditional  mean 
scalar  diffusion  rate,  and  the  mixture  fraction  probability  density  function  (PDF)  were 
determined  from  the  statistically  stationary  simulation  data. 

The  mean  scalar  dissipation  rate  (Nv  =  (Z?(V£)2  |  £(z, t)  =  77)),  conditional  upon  the 
value  (77)  of  the  local  instantaneous  value  of  mixture  fraction  (£),  was  found  to  be  inde¬ 
pendent  of  mixture  fraction  as  reported  by  Overholt  and  Pope[7],  and  Pope  and  Ching[8]. 
Thus  the  significant  simplification, 

Nv  =  (N\r])  =  (N)  ,  (5) 

was  employed  in  the  modelling  described  below. 

The  mixture  fraction  PDF  (Pv)  for  the  whole  domain  was  found  to  agree  closely  with 
a  superposition  of  Gaussian  PDFs  with  equal  variance,  but  with  mean  values  varying 
with  position  in  the  x  direction  according  to  the  imposed  mean  mixture  fraction  gradient. 
Given  values  of  mean  mixture  fraction  on  the  x  boundaries  of  £0  and  £1,  and  a  spatially 
uniform  mixture  fraction  variance  a2,  the  mixture  fraction  PDF  can  be  written  as, 
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where  <f>(. . .)  is  the  integral  of  the  Gaussian  distribution  between  the  given  argument  and 
infinity. 

The  conditional  mean  scalar  diffusion  rate  (Mv  =  (PV2£  |  £(z,  t)  =  rj))  is  yielded  by, 


■ 


which  holds  for  homogenous  mixing  conditions.  The  simplification  afforded  by  Eqn  5,  thus 
gives  Mv  as, 

„  _  W  ap "  IS) 

M "  -  nT  ay  •  (8) 

The  conditional  diffusion  profile  given  by  the  equation  above  was  found  to  closely  agree 
with  the  simulation  data. 


Figure  1:  Profiles  for  mixture  fraction  PDF  (O),  conditional  mean  scalar  dissipation  rate 
(Q),  and  conditional  mean  diffusion  rate  (A),  expected  and  observed  in  simulation  data. 

The  forms  of  the  mixture  fraction  PDF  and  conditional  diffusion  rate  profiles,  given 
by  the  above  equations,  are  plotted  in  Fig.  1.  The  mixture  fraction  PDF  is  somewhat 
broadened  around  the  peak  compared  to  any  single  Gaussian  profile  with  the  same  mean 
and  variance.  As  a  result,  the  conditional  diffusion  rate  profile  has  a  mean  slope  which 
is  about  35%  shallower  than  that  of  the  purely  linear  profile  which  results  from  a  single 
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Gaussian  PDF  with  an  equivalent  mean  and  variance,  and  equivalent  conditional  mean 
dissipation  rate. 

All  of  the  simulations  exhibited  profiles  of  the  nature  depicted  in  the  figure,  since  they 
all  shared  exactly  the  same  mixing  characteristics.  The  spatially  uniform  value  of  variance 
of  mixture  fraction  was  equal  to  ~  0.34  in  each  case,  with  mean  mixture  fractions  on  the 
x  boundaries  of  zero  and  unity,  and  an  overall  domain-averaged  mean  mixture  fraction  of 
one  half. 

Notice  from  Fig.  1  that  the  quantity  £,  refered  to  here  as  mixture  fraction ,  is  normalised 
to  have  a  unit  value  on  the  rich  boundary  of  the  domain  and  zero  value  on  the  lean 
boundary,  but  this  does  not  denote  pure  fuel  and  oxidizer  states.  It  is  possible,  and  indeed 
required  in  the  current  simulation  configuration,  to  have  values  of  £  which  are  greater  than 
unity  and  less  than  zero.  Pure  fuel  and  oxidizer  states  can  be  considered  to  exist  only  in 
the  limits  of  £  — »  ±oo. 


3  Modelling  Method 

Application  of  the  CMC  model  for  gas  phase  turbulent  combustion  to  the  modelling  of 
condensed-phase  fields  requires  some  consideration  of  the  differences  in  behaviour  between 
the  two  phases.  Firstly,  it  was  assumed  that  the  particles  were,  in  the  mean,  small  enough 
to  follow  the  smallest  scales  of  motion.  This  was  true  in  the  inertia-less  particle  simulations, 
but  was  somewhat  less  justified  for  the  inertial  (particularly  droplet)  simulations. 

Secondly,  it  is  clear  that  the  particle  phase  does  not  diffuse  appreciably  on  a  molecular 
level,  and  thus  the  particle  phase  and  the  gas  phase-based  mixture  fraction  are  transported 
in  a  completely  different  manner  at  this  level.  The  effect  of  this  differential  diffusion  must 
therefore  be  incorporated  into  the  CMC  model. 

Further,  at  the  smallest  scales,  a  distinct  difference  between  the  distribution  of  gaseous 
species  and  condensed  phase  species  is  expected  (see  Klimenko[l]).  Significant  local  struc¬ 
ture  can  exist  in  the  particle  field  at  scales  where  gas  phase  fields  have  been  completely 
smeared  by  molecular  diffusion.  The  existence  of  such  local  structure  implies  that  instan¬ 
taneous  local  deviations  of  particle  mass  fractions  from  means,  conditional  upon  mixture 
fraction,  can  be  large. 

In  the  present  study,  particle  effects  upon  the  gas  phase  were  neglected,  leaving  the 
effect  of  gas  phase  fluctuations  on  particle  evolution  to  be  examined  in  isolation.  In  reality, 
however,  reciprocal  interactions  between  gaseous  and  condensed-phase  species  is  likely  to 
be  very  important.  In  the  case  of  droplet  evaporation,  it  is  the  addition  of  vaporised  fuel 
to  the  continuous  phase  which  largely  determines  the  way  in  which  combustion  proceeds. 
In  the  case  of  soot  particles,  the  removal  of  sensible  enthalpy  from  the  continuous  phase 
occurs  as  a  result  of  proximity  to  soot,  which  can  strongly  affect  localized  combustion 
dynamics.  In  both  of  these  instances,  large  conditional  variance  in  particle  properties 
could  lead  to  similarly  elevated  levels  of  conditional  variance  in  sympathetic  gas  phase 
species  and  thereby  increase  the  difficulty  of  the  chemical  closure  problem  in  the  gas 
phase.  This  increase  in  conditional  variance  for  gas  phase  species  can  be  viewed  as  arising 
from  the  different  gas  phase  behaviour,  at  the  same  mixture  fraction,  which  will  result 
depending  on  whether  a  parcel  of  gas  is  adjacent  to  a  particle  or  not. 
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It  is  for  this  reason  that  it  is  important  to  be  able  to  predict  the  level  of  conditional 
variance  of  particle  mass  fraction.  The  CMC  model  proposed  in  this  study  made  use 
of  a  particle-specific  differential  diffusion  variant  of  the  second-order  conditional  moment 
closure  proposed  by  Li  and  Bilger[9]. 


3.1  Model  Derivation 


In  the  following,  a  CMC  model  for  particle  reactions  in  turbulent  nonpremixed  com¬ 
bustion  is  derived  for  general  flow  and  mixing  conditions.  The  spatial  and  temporal 
simplifications,  afforded  by  the  simulation  conditions  of  the  present  study,  are  introduced 
in  Section  3.2. 


The  local  instantaneous  equations  for  mixture  fraction  and  a  particle  mass  fraction 
continuum  form  the  basis  for  the  derivation  of  the  model.  The  equation  for  mixture 
fraction  £  can  be  written  as, 


n 

dt+Uldxi  vdx?  0 


(9) 


where  the  flow  field  is  assumed  to  be  incompressible,  U{  is  the  component  of  carrier  fluid 
velocity  in  the  2-th  direction,  and  V  is  the  diffusivity  of  mixture  fraction  under  a  Fickian 
transport  assumption. 


by, 


The  equation  for  local  instantaneous  mass  fraction  of  the  particle  continuum  is  given 


dY  8Y  d2Y  _  . 
dt  Vl  dxi  v  dx2-  U 


(10) 


where  w  is  the  local  reaction  rate  of  the  particle  continuum,  Dy  is  the  effective  diffusivity 
of  this  continuum,  and  V{  is  the  velocity  component  in  the  i-th  direction,  of  the  field  of 
particles.  The  particle  diffusivity  Vy  is  expected  to  be  very  much  smaller  than  the  mix¬ 
ture  fraction  diffusivity,  however,  the  term  is  retained  in  the  derivation  for  completeness. 
Multiplying  Eqn  10  by  the  particle  mass  fraction  yields, 


dY 2  BY2  w  fl2Y 

+  Vi— - 2D„Y  -rr-r  =  2Yw 


dt 


dx. 


dx2 


(11) 


which  is  used  in  the  derivation  of  the  equation  for  conditional  variance  of  particle  mass 
fraction. 


The  fine-grain  probability  density  function  (see  Pope[7])  for  mixture  fraction  is  defined 
as, 

ip{x,t,rj)  =  5(^{x,t)-ri)  ,  (12) 

and  has  an  expectation  over  the  statistical  ensemble  equal  to  the  mixture  fraction  PDF 
(Pv)  as  given  by, 

P-n  (z>  t,  V)  =  W  •  (13) 

The  local  instantaneous  equation  for  the  fine-grained  PDF  can  be  derived  from  the  differ¬ 
ential  properties  of  if)  and  Eqn  9,  so  that, 

W  +  Wi  ^  +  Wn  (®,6v2f)  =  0  '  <14) 
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Klimenko  and  Bilger[10]  derived  a  conditional  mean  equation  for  the  differential  diffu¬ 
sion  of  gaseous  species  in  turbulent  combustion  through  the  combination  of  equations  11 
and  14  given  above.  Their  methodology  is  largely  followed  here,  except  that  a  conditional 
variance  equation  is  also  derived  and  some  significant  differences  in  closure  assumptions  are 
made  in  the  final  stages.  In  the  section  immediately  following,  derivation  of  the  equation 
for  the  product  of  the  square  of  particle  mass  fraction  and  the  fine  grain  PDF  is  described, 
as  a  prerequisite  to  the  derivation  of  the  conditional  variance  equation  for  particle  mass 
fraction.  The  derivation  of  the  corresponding  equation  for  conditional  mean  mass  fraction 
is  analogous,  except  where  noted  otherwise. 


3.1.1  Derivation  of  ipY2  equation 


The  derivation  proceeds  by  adding  Eqn  11,  multiplied  by  ip,  to  Eqn  14,  multiplied  by 
Y2.  The  resultant  equation  is  given  by, 

n  c\  O'V/2 

Ft  (^y2)  +  dxl  (Ui'^y2)  =  (Ul  ~  Vi)  ^~FF  +  2^YCj  +  Av~At  >  (15) 

where 

.  a(  =  , 

are  convenient  groupings  of  the  V  and  Vy  diffusive  terms.  The  mixture  fraction  diffusivity 
term,  first  expanded  by  using  the  inhomogeneous-flow  form  of  Eqn  7,  is 

A£  =  y2(l?(mv?)VI,^«’))  ■  <16) 

The  latter  right-hand  side  term  appearing  in  Eqn  16,  also  appears  in  the  full  expansion  of 
{Y2ip).  Substitution  of  the  expansion  into  Eqn  16  yields  a  final  expression  for  A$, 


The  particle  diffusivity  term  Ay,  as  given  above  in  connection  with  Eqn  15,  appears 
in  the  full  expansion  of  VyY  §|-  (*&)■  Expressing  Ay  as  the  subject  of  this  expansion, 


yields 


"  di 7 


d£  dY2\  ^  „  8  (  ,  dY\ 
dxj  dxj  )  +WyYdxj  ydxj)  ' 


(18) 


Rearrangement  of  the  latter  right-hand  side  term  of  the  above  equation  via  the  chain  rule 
gives  Ay  in  final  form, 


A,, 


d_ 

drj 


Vylp 


dt e  dY2' 

dxj  dxj 


d_ 

1  dxj 


+  vy^~  ^ 


dY2' 

dx. 


~  2  Vylp 


dY_ 

dxi 


The  substitution  of  Eqns  17  and  19  into  Eqn  15  yields, 


+  ^w+Rn 


(19) 


(20) 
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where  the  flux  of  squared  particle  mass  in  mixture  fraction  space  is  given  by, 

Fy2  =  {Vt  +  V)i,K^.^(Y^V{V(?)  ,  (21) 

and  the  collected  residual  terms  can  be  written  as, 

HY2  f)  (  FfY 2  \  ft2  /  \ 

Rn  „  („  -  +  (V,  -  V)  Wj  ^ J  +  Bjg  (y  V)  •  (22) 

3.1.2  Averaging  and  closure  of  xpY  equation 

It  is  illustrative  to  average  and  close  the  simpler  equation  for  conditional  mean  particle 
mass  fraction  before  employing  the  same  approach  on  the  conditional  variance  equation. 
The  conditional  mean  equation  is  required  before  further  derivation  of  the  conditional 
variance  equation  can  proceed.  In  this  section,  some  of  the  key  closure  assumptions  of 
Klimenko  and  Bilger[10]  have  been  relaxed  out  of  neccessity  for  the  modelling  of  particles. 

A  set  of  equations,  similar  to  Eqns  20-22,  apply  for  the  product  of  the  fine-grain  PDF 
and  particle  mass  fraction,  and  can  be  written  as, 

|Wy)  +  |-(^y)  =  ^+^  +  Kr,  (23) 

Fy  =  (V,  +  V)  {y#D  (VQ2)  .  (24) 

and 

Ry  =  (*  -  *)  ^  +  (B,  -  V)  A  (*Q  +  or*)  ■  (25) 

Averaging  Eqn  23  over  mixture  fraction  space  yields  the  PDF-product  form  of  the  condi¬ 
tional  moment  closure  (CMC)  equation  for  conditional  mean  particle  mass  fraction, 

Jt +  £ {{uiY  1  ’I>  F")  =  <" 1  n)  p,  +  ^  >  (26) 


where  Qv  =  (Y  \  rj)  =  (Yip)  /Pv,  and  Fq  and  Rq  are  the  conditionally  averaged  flux  and 
residual  terms.  The  averaged  residual  term  can  be  written  as  follows, 


and  contains  terms  which  are  assumed  to  be  small  for  high  Reynolds  number  flow  with 
only  small  particle  slip  velocities.  The  conditional  mean  flux  term,  equal  to 


J%-(P,  +  D)P,(^|,)-A(l.P,(y(VQ>|,»  .  (28) 

is  approximated  by  Klimenko  and  Bilger[10],  after  assuming  that  the  gradient  of  Y  is  well 
correlated  with  the  gradient  of  mixture  fraction,  and  that  there  is  only  a  weak  correlation 
between  Y  and  (V£)2.  In  the  following,  the  latter  assumption  is  relaxed,  so  that 

S3«cP,JV„^-0,^(R,JV,)-^(F„C,)  ,  (29) 
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where  Nv  =  {V  (V£)2  |  rjj  is  the  scalar  dissipation  rate,  and  C n  is  the  covariance  between 
scalar  dissipation  rate  and  particle  mass  fraction.  The  symbol  c  denotes  a  mixing  mode 
variable  which  is  equal  to  unity  in  the  case  of  mixing  in  the  absence  of  differential  diffusion, 
but  is  not  known  for  the  case  of  differentially  diffusing  scalars  (Klimenko  &:  Bilger  [10]). 
This  variable  is  discussed  in  greater  detail  in  Section  3.2. 

The  final  form  of  the  equation  for  conditional  mean  particle  mass  fraction  is  obtained 
by  substituting  Eqn  29  into  Eqn  26,  followed  by  subtraction  of  the  mixture  fraction  PDF 
equation  (averaged  form  of  Eqn  14),  and  normalisation  by  the  PDF,  to  give, 


dQr, 

dt 


+  (ui  |  v) 


dQv 

dxi 


=  (W  I  V)  +  cNr, 


d'Qr, 

dr/ 2 


+  (c  —  \)Mn 


dQv 

dr) 


1  d 2 


ORA)  • 


(30) 


In  the  above  expression,  Rq  has  been  neglected,  as  has  the  conditional  correlation  between 
velocity  fluctuations  and  particle  mass  fraction.  As  in  Eqn  7,  the  symbol  Mv  denotes  the 
conditional  mean  diffusion  rate. 


3.1.3  Averaging  and  closure  of  ip Y 2  equation 


The  equation  for  the  conditional  variance  of  particle  mass  fraction,  defined  a s  <?2  = 
( y 2  |  rj)  =  (y ’2ip)  /P^  -  Q 2,  can  be  derived  from  Eqn  20  through  several  stages.  Averaging 
of  Eqn  20  yields  a  mean-square-PDF  product  equation, 

J)  («,  {Y2  b»  +  ((«.'  V2  I  v)  P„)  =  2  (Yu  |  v)  P,  +  j*  -  p,er  +  R,  .  (31) 

where  Fq  denotes  the  flux  of  mean  square  mass  in  mixture  fraction  space,  Rq  is  the  average 
of  the  y2  residual  term  of  Eqn  22, 

r,  =  P,  («,  -  *)  I  v)  +  fl>,  -  V)  £-  (p,  (gf  I  *))  +  (p,  (r2 1 i»  . 

(32) 

and  ey  is  the  conditional  mean  square  mass  dissipation  rate,  given  by 

(y  =  2Vy((VY)2  |  t?)  .  (33) 

Following  the  closure  argument  for  flux  of  conditional  mean  mass  in  mixture  fraction  space 
(Eqn  29),  the  flux  of  the  conditional  mean  square  particle  mass  fraction, 

F„  =  (®,  +  V)  P,  (£§£  I  ,)  -  £  (P,  (Y2V  ( V?)’  I  ,»  ,  (34) 


is  approximated  by  the  following,  where  the  mixing  mode  variable  c  has  a  value  equal  to 
that  employed  in  Eqn  29, 


dg2 


c  -  q2-^  ( PVN„ )  +  2c  Pvcn^  -  2 (PnCv) 


dQr, 


dr) 


8_ 

1  drj 


PcPrjNrj  Qlflr,  (PriNv)  ■ 


dr) 


(35) 
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Following  the  substitution  of  Eqn  35  into  Eqn  31,  the  decomposition  of  Y  into  Qn  +y, 
and  the  subtraction  of  the  equation  for  the  conditional  mean  mass  fraction,  the  equa¬ 
tion  for  the  product  of  conditional  variance  and  mixture  fraction  PDF  is  derived.  The 
decomposition  of  this  equation  via  the  chain  rule,  and  subtraction  of  the  PDF  equation 
multiplied  by  q2  yields  the  final  form  of  the  conditional  variance  equation, 


+  (“*’  I  ^  =  2  ^  I  ^  +  cNr,~dr^  +  (c  “  !)  ’  (36) 


dx; 


where  Sq  denotes  the  collected  conditional  source  terms  for  particle  mass  fraction  variance, 


s  ^e~w +hc-1]r,  +*(«-£) «.  (w)'' 


and  e„  is  the  dissipation  rate  of  conditional  variance, 


=  2Vy((Vy)2\ri)  . 


(37) 


(38) 


Note  that  the  conditional  variance  equation,  in  the  form  given  above,  neglects  conditional 
covariance  between  velocity  and  mass  fraction,  and  all  of  the  terms  associated  with  the 
difference  between  Rq  and  2QvRq  which  appear  in  the  derivation. 

The  final  forms  of  the  equations  for  the  conditional  mean  and  variance  of  particle  mass 
fraction  (Eqns  30  &  36)  contain  two  unclosed  quantities,  namely  the  mixing  mode  variable 
c,  which  parameterizes  the  level  of  differential  diffusion,  and  the  covariance  between  scalar 
dissipation  rate  and  particle  mass  fraction  Cv.  A  discussion  of  these  two  quantities  is 
provided  in  the  following  sections. 


3.2  Solution  of  Model  Equations 

Equations  30  &  36,  of  the  preceding  section,  were  solved  to  match  the  simulated  condi¬ 
tions  of  spatially  homogeneous,  steady  conditional  statistics.  The  statistically  stationary, 
spatially  degenerate  form  of  the  conditional  mean  particle  mass  fraction  profile  is  given 

by, 

B20  BO  1  B 2 

0={u\v)  +  eNv-^  +  (c-l)Mvj^--—2(PriCv)  +  lQ  ,  (39) 

where  Iq  is  the  conditional  mean  particle  mass  injection  rate  which  results  from  the  injec¬ 
tion  process  described  in  Section  2.1.  The  conditional  variance  of  particle  mass  fraction 
profile  was  found  from, 

BP'q2  Bq^ 

0  =  2(yO'\r1)  +  cNv-^-  +  (c-l)M7]-^  +  Sq  +  Iq  ,  (40) 

where  Iq  is  the  conditional  variance  source  due  to  particle  injection.  Note  that  the  effective 
diffusivity  of  the  particle  field  (Vy)  was  assumed  to  be  identically  zero,  thus  the  absence  of 
the  eq  term,  given  in  Eqn  36,  from  the  above  equation.  The  conditional  particle  injection 
source  terms  included  in  the  above  (Iq  &  Ig )  are  given  by, 

Iq  =  ~~  (1  —  Qn)  (41) 

•pop 
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r,  =  {(i  -  Q,)2  -  it)  (42) 

T pop 

where  Pj  is  the  PDF  of  particle  injection  mixture  fraction,  and  Tpop  is  the  time  required  for 
the  number  of  injections  to  equal  the  size  of  the  particle  population.  In  the  above,  it  was 
assumed  that  particles  were  injected  to  replace  those  which  leave  the  system,  and  all  freshly 
injected  particles  had  unit  mass,  as  described  in  Section  2.1.  The  form  of  the  injection 
PDF  consisted  of  the  normalised  sum  of  two  Gaussian  distributions  with  variances  equal 
to  the  domain  variance,  and  mean  mixture  fraction  values  of  zero  and  unity. 

Chemical  source  terms  in  Eqns  39  &  40  were  determined  using  instantaneous  local 
temperature  (taken  from  an  invariant  chemical  equilibrium  profile),  and  the  conditional 
mean  and  variance  of  particle  mass  fraction  in  truncated  Taylor  series  expansions  of  the 
instantaneous  reaction  rate  expressions. 

The  mixing  mode  variable  c,  included  in  the  CMC  equations  above,  provides  a  means 
of  accounting  for  differential  diffusion  effects.  In  the  limit  where  mixing  is  overwhelmingly 
due  to  turbulent  stirring,  the  mixing  mode  variable  tends  to  unity.  In  cases  where  a  larger 
proportion  of  the  mixing  is  due  to  molecular  mass  transfer,  the  mixing  mode  variable  will 
vary  substantially  from  unity.  The  greatest  allowable  deviation  of  c  from  unity  is  given 
by  the  ratio  of  the  molecular  diffusivity  of  the  species  in  question,  to  that  of  the  mixture 
fraction. 

From  a  modellers  perspective,  it  is  highly  desirable  that  c  be  independent  of  Qv,  q*.  and 
mixture  fraction.  In  that  eventuality,  one  could  hypothesize  that  the  degree  of  departure, 
from  mixing  that  is  dominated  by  strong  turbulent  stirring,  might  be  expressed  as, 

C=l+(|-l)  fmix  ,  (43) 

where  /TOjr  is  a  function  which  varies  between  zero  and  unity  according  to  some  dependence 
on  global  properties  of  the  flow  field.  One  possibility  for  fmix  is  some  power  of  the 
Kolmogorov  scalar  scale  £*,  appropriately  normalised  by  the  maximum  globally  realizable 
range  in  mixture  fraction.  The  Kolmogorov  scalar  scale  is  a  measure  of  the  characteristic 
size  of  scalar  fluctuations  at  dissipative  scales.  For  mixture  fraction,  this  scale  varies 
between  zero  and  unity  for  high  and  low  intensity  turbulent  mixing  respectively. 

As  this  study  consists  of  data  from  a  single  set  of  mixing  conditions,  it  is  not  possible 
to  validate  any  model  for  variation  in  c  with  mixing  conditions.  It  is,  however,  possible  to 
examine  whether  a  single  value  of  c  is  appropriate  independently  of  local  values  of  Qv,  q%, 
and  mixture  fraction.  This  issue  is  examined  in  the  next  section. 

It  should  be  noted  that  Kronenburg  and  Bilger[ll]  successfully  account  for  differ¬ 
ential  diffusion  in  the  CMC  model  through  the  use  of  an  alternate  treatment  to  that 
proposed  above.  Their  methodology  has  much  to  recommend  it,  but  when  applied  to 
multi-component  reacting  systems,  it  requires  an  additional  conditional  moment  equation 
to  be  solved  for  each  species.  The  method  presented  here  seeks  to  avoid  that  added  cost. 

It  was  found  that  the  profile  shape  and  magnitude  of  the  conditional  covariance  between 
particle  mass  fraction  and  scalar  dissipation  rate  ( Cv )  was  critical  to  the  transport  of 
particles  in  mixture  fraction  space.  While  the  form  of  this  profile  is,  in  general,  not 
known,  its  behaviour  at  mixture  fraction  bounds  is  prescribed  by  the  integration  of  Eqn 
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□  Covar-Pred 


Figure  2:  Modelled  and  observed  conditional  covariance  between  mass  fraction  and  scalar 
dissipation  rate  for  case  si.  Symbols  denote  the  following  profiles  :  □  -  prediction  of 
covariance,  and  A  -  simulated  covariance. 
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26,  with  respect  to  mixture  fraction,  between  these  bounds.  This  integration  should  yield 
the  unconditional  mean  mixture  fraction  equation,  which  implies  that  the  bounding  values 
of  the  mixture  fraction  flux  term  (Eqn  29)  are  zero.  It  follows  that,  the  slope  of  the  product 
must  be  zero  at  the  bounding  mixture  fractions. 

Modelling  of  the  Cn  term  was  achieved  using  a  simple  heuristic  approach,  to  best  match 
the  observed  simulation  profiles  for  Cn.  This  model  took  the  form, 

Crj  =  —  llQQrj  +  Ctcor^  -^7)  >  (44) 

where  aq,  aQ  and  acor  were  constants  determined  from  the  goodness  of  fit  with  the  Cv 
profile  observed  from  the  corresponding  simulation.  The  above  fitting  technique  is  similar 
to  that  employed  by  Li  and  Bilger[9],  except  for  the  inclusion  of  the  conditional  mean  value 
(Qn),  and  the  zero-boundary  correction  term  acor.  It  was  found  that  it  was  necessary  to 
include  these  terms  to  match  the  observed  simulation  data.  Due  to  the  case-specific  nature 
of  the  covariance  profiles,  they  are  discussed  in  turn  in  the  following  section. 


4  Results 


In  the  following,  the  results  from  the  modelling  of  soot  particles  and  droplets  are 
discussed  separately. 


4.1  Soot  particle  results 

Examining  the  conditional  covariance  between  scalar  dissipation  rate  and  particle  mass 
fraction  {Cv),  the  simulation  data  revealed  the  characteristic  form  of  the  covariance  profile 
for  soot  as  possessing  a  large  negative  spike  near  the  mean  mixture  fraction.  The  location  of 
this  spike  in  mixture  fraction  space  corresponded  with  the  location  of  the  peak  conditional 
mean  mass  fraction.  Through  careful  selection  of  the  coefficients  in  the  Cn  model  equation 
(Eqn  44),  it  was  possible  to  produce  the  similar  profile  as  wras  observed  in  the  simulation 
data,  while  producing  good  predictions  for  Qv  and  q 

The  profile  match  derived  for  conditional  covariance  in  soot  case  si  can  be  seen  in  Fig¬ 
ure  2.  Note  that  the  observed  simulation  profile  has  been  slightly  smoothed,  and  averaged 
about  the  symmetry  point  at  the  mean  mixture  fraction,  to  facilitate  comparison  with  the 
smooth  symmetric  model  profile.  It  is  evident  that  the  model  profile  agrees  reasonably  well 
with  the  observed  profile,  particularly  in  the  vicinity  of  the  mean  mixture  fraction.  The 
model  profile  does  not  match  the  observed  maximum  value  of  the  covariance  at  mixture 
fractions  slightly  leaner  and  richer  than  the  mean  value.  Unfortunately,  it  is  the  first  and 
second  derivatives  of  the  conditional  covariance  which  appear  in  the  model  equations.  The 
value  of  the  slope  and  curvature  of  the  model  C n  profile  are  believed  to  reflect  the  general 
behaviour  of  the  simulation  data,  however  it  is  difficult  to  make  any  stronger  statement 
due  to  the  limited  data  available  for  statistical  analysis.  Future  simulations  with  larger 
particle  populations  will  be  used  to  better  understand  this  behaviour. 

The  characteristic  effect  on  the  conditional  mean  particle  mass  fraction  of  varying  the 
mixing  mode  variable,  c,  can  be  seen  in  Figure  3  for  case  si,  along  with  the  conditional 


14 


DSTO-RR-0 142 


A  c=005 


Figure  3:  Predicted  conditional  mean  mass  fraction  profiles  for  case  si  with  arbitrarily 
varied  levels  of  differential  diffusion.  Symbols  denote  the  following  profiles  :  A  -  prediction 
with  c  =  0.05,  o  -  prediction  with  c  =  1.00,  and  o  -  chemical  equilibrium  with  particle 
injection. 
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Mixture  Fraction 


Figure  4:  Predicted  conditional  root  mean  square  deviation  mass  fraction  profiles  for  case 
si  with  arbitrarily  varied  levels  of  differential  diffusion.  Symbols  denote  the  following 
profiles  :  A  -  prediction  with  c  =  0.05,  o  -  prediction  with  c  —  0.50,  and  o  -  prediction 
with  c  =  1.00. 
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mean  profile  which  results  from  chemical  reactions  and  particle  injection  alone.  It  can  be 
seen  that  the  arbitrary  increase  of  c  causes  the  predicted  mean  mass  fraction  profile  to 
simultaneously  increase  in  peak  value  and  decrease  in  minimum  value.  Further,  increases 
in  c  tend  to  decrease  the  value  of  conditional  mean  mass  fraction  at  very  rich  and  very  lean 
mixture  fractions.  The  governing  conditional  mean  particle  mass  fraction  equation  (Eqn 
39),  has  the  property  that  increased  values  of  c  tend  to  favor  a  mode  of  mass  transfer  in 
mixture  fraction  space  which  is  essentially  diffusive  in  nature.  Lower  values  of  c  favor  the 
mode  associated  with  the  drift  of  particle  mass  in  mixture  fraction  space  which  results  from 
gas  phase  species  diffusing  at  a  more  rapid  rate  in  physical  space  than  the  particles  they 
surround.  Thus  c  can  be  thought  of  as  a  measure  of  to  what  degree  particle  mass  diffuses 
against  a  mixture  fraction  coordinate,  instead  of  undergoing  an  apparent  convection-like 
process  in  mixture  fraction  space  due  to  the  differential  diffusion  of  the  coordinate  and 
mass  fraction  fields  in  physical  space.  The  tendency  towards  smoother  conditional  mean 
profiles,  in  Fig.  3,  with  increased  c  values  is  indicative  of  this  trend. 

Variation  in  the  predicted  conditional  mass  fraction  variance  profiles  with  arbitrary 
variation  in  the  mixing  mode  variable  can  be  seen  in  Figure  4  for  case  si.  It  is  apparent 
that  increasing  the  value  of  c  leads  to  an  increased  peak  value  in  the  conditional  mean 
variance  profiles  and  an  increased  value  of  variance  at  very  rich  and  lean  mixture  fractions. 
The  elevated  level  of  conditional  variance  in  the  high  c  case  is  due  to  the  source  contribution 
from  the  last  term  of  Sq  (see  Eqn  37). 

A  low  value  of  c  results  in  a  dip  in  variance  near  the  mean  mixture  fraction,  where  the 
model  equations  do  not  predict  any  variance  production.  This  dip  is  not  present  in  higher 
c  cases  where  diffusive  transport  in  mixture  fraction  space  is  strong  enough  to  smooth  out 
sharp  changes  in  gradient. 

The  result  of  matching  the  conditional  covariance  profile  and  minimizing  the  mixing 
mode  variable  to  c  =  0.05  (minimum  numerically  stable  value)  provided  the  predicted 
conditional  mean  and  variance  profiles  plotted  in  Figure  5. 

Employing  the  same  value  of  c  and  the  same  constants  in  the  model  for  the  mass- 
dissipation  covariance  Cv  from  Eqn  44,  provided  accurate  predictions  for  the  conditional 
mean  and  variance  of  particle  mass  in  soot  case  s2i.  These  predictions  are  plotted  against 
simulation  data  in  Figure  6. 

It  is  evident  that  there  is  reasonable  qualitative  agreement  between  the  predicted  and 
observed  profiles  in  Figs  5  &  6.  The  CMC  model  captures  the  essential  changes  in  profile 
shape  in  the  more  reactive  case  (Fig.  5),  but  appears  to  incorrectly  predict  rich  and  lean 
side  behaviour  to  greater  extent  in  the  less  reactive  case.  Good  quantitative  agreement 
is  found  near  the  mean  mixture  fraction  in  both  cases,  however  the  predicted  profiles 
deviate  substantially  from  the  observed  data  at  mixture  fractions  away  from  the  mean. 
This  discrepancy  may  have  been  due  to  the  poor  modelling  of  the  conditional  covariance 
term  ( Cn )  at  those  same  mixture  fractions.  It  is  difficult  to  be  sure  that  the  model  fit 
to  Cv  provided  by  Eqn  44,  provides  accurate  Cv  gradients  and  curvatures  as  is  required 
in  the  CMC  equations.  It  is  clear  that  the  modelling  of  the  mass-dissipation  covariance 
requires  further  attention. 
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Figure  5:  Modelled  and  observed  conditional  mean  and  variance  of  particle  mass  fraction 
for  case  si.  Symbols  denote  the  following  profiles:  A  -  Qv  prediction,  o  -  simulated  Qn,  o 
-  qv  prediction,  and  x  -  simulated  qn. 
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Figure  6:  Modelled  and  observed  conditional  mean  and  variance  of  particle  mass  fraction 
for  case  s2i.  Symbols  denote  the  following  profiles:  A  -  Qv  prediction,  o  -  simulated  Qv, 
o  -  qv  prediction,  and  x  -  simulated  q v. 
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4.1.1  Effect  of  particle  slip  velocity 

Soot  cases  s2  and  s2i  were  simulated  to  examine  the  effect  of  slip  velocities  between 
particles  and  their  carrier  fluid.  The  former  simulated  case  took  account  of  particle  inertia, 
and  the  difference  in  local  velocities  which  can  arise  from  the  kinetic  lag  between  fluid  flow 
variations  and  particle  response  to  those  variations.  The  latter  simulation  case  ( s2i )  did 
not  account  for  particle  inertia  and  thus  no  slip  velocities  were  present. 

The  differences  in  the  conditional  mean  and  variance  particle  mass  fraction  profiles 
between  these  two  simulation  cases  can  be  seen  in  Figure  7.  It  is  evident  that  the  intro¬ 
duction  of  particle  inertia,  all  else  being  equal,  causes  the  conditional  mean  particle  field 
to  behave  in  a  more  diffusive  manner  in  mixture  fraction  space.  Simultaneously,  a  rise  in 
conditional  variance  results,  particularly  at  mixture  fractions  away  from  the  mean  mix¬ 
ture  fraction.  Although  not  plotted,  there  appeared  to  be  little  change  in  the  Cv  profiles 
between  the  simulation  pair,  indicating  that  the  vee-shaped  form  of  the  covariance  profile 
is  not  a  result  of  inertial  particle  transport. 

Conditional  source  terms  related  to  particle  slip  velocity  appear  in  the  neglected  resid¬ 
ual  groupings  of  Eqns  39  k.  40.  The  conditional  mean  equation  contains  the  slip  trans¬ 
port  term  ((u,-  —  v,)  |  rf)  ^  {Qv),  whereas  the  variance  equation  contains  a  transport  term 

(( u\  -  v')  |  rj )  and  the  source  term  ((u'-  -  v[)y  |  t?)  ~  ( Qv ).  The  difference  in  be¬ 

haviour  observed  in  the  above  simulation  cases  corresponds  with  the  expected  effect  of 
the  neglected  slip  velocity  terms  in  the  CMC  equations.  It  is  clear  that  the  introduction 
of  larger  inertial  particles  to  the  model  problem  gives  rise  to  a  significant  increase  in  the 
difficulty  of  modelling  particle  behaviour. 


4.2  Droplet  particle  results 

Unlike  the  soot  particles,  the  droplet  particle  simulations  exhibited  very  weak  condi¬ 
tional  covariance  between  scalar  dissipation  rate  and  particle  mass  fraction.  This  distinc¬ 
tion  between  the  particle  types  suggests  that  the  form  of  the  covariance  (Cj,)  profile  is 
linked  to  the  types  of  reactions  which  act  upon  the  particles.  It  would  seem  that  strong 
reactions,  as  in  the  soot  cases,  give  rise  to  sharp  changes  in  particle  properties  in  mix¬ 
ture  fraction  space  which  in  turn  impact  strongly  upon  the  transport  statistics  between 
adjacent  mixture  fractions.  The  weak  reactions  of  the  droplet  cases  do  not  impart  sharp 
changes  in  value  to  the  conditional  statistical  profiles,  and  thus  do  not  seem  to  cause  strong 
covariance  between  mass  fraction  and  scalar  dissipation  rate.  Further  investigation  of  this 
behaviour  is  clearly  warranted. 

Due  to  the  observed  behaviour  in  the  droplet  cases,  the  covariance  profile  was  modelled 
as  being  zero  for  all  mixture  fractions.  As  with  the  soot  cases,  the  value  of  the  mixing 
mode  parameter,  which  gave  best  predicted  agreement  with  the  observed  mass  fraction 
profiles,  was  very  low  (c  =  0.1).  The  degree  of  agreement  between  the  observed  and 
predicted  profiles  for  the  inertia-less  droplet  case  (dli)  can  be  inferred  from  Figure  8.  It 
is  clear  that  the  CMC  model  predicts  the  conditional  mean  and  variance  of  droplet  mass 
fraction  with  a  good  degree  of  accuracy. 
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Figure  7:  Observed  conditional  mean  and  variance  of  particle  mass  fraction  for  inertial  (s2) 
and  inertia-less  (s2i J  soot  cases.  Symbols  denote  the  following  profiles:  A  -  Q v  inertial,  o 
-  Qv  inertia-less,  o  -  qv  inertial,  and  X  -  qr.  inertia-less. 
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Figure  8:  Modelled  and  observed  conditional  mean  and  variance  of  particle  mass  fraction 
for  case  dli.  Symbols  denote  the  following  profiles:  A  -  Qv  prediction,  o  -  simulated  Qv, 
o  -  qv  prediction ,  and  x  -  simulated  qn. 
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4.2.1  Effect  of  particle  slip  velocity 

While  it  is  clear  that  the  CMC  model  can  capture  the  behaviour  of  very  small  fuel 
droplets  in  the  terminal  stage  of  evaporation,  the  modelling  of  particles  with  significant 
inertia  is  more  problematic. 

In  Figure  9,  the  conditional  statistics  from  the  simulation  pair  of  dl  and  dli  are  com¬ 
pared.  It  is  evident  that,  as  with  the  soot  comparison  of  Fig.  7,  the  larger  inertial  particles 
exhibit  a  tendency  towards  conditional  mean  profiles  which  have  a  lower  level  of  curvature, 
as  though  smoothed  by  enhanced  diffusive  transport.  The  conditional  variance  profile  ex¬ 
hibits  virtually  no  change  in  slope  at  all  in  the  case  of  the  inertial  particles.  The  existence 
of  slip  velocities,  larger  than  in  the  inertial  soot  case,  causes  a  great  deal  of  conditional 
variance  generation  to  occur  at  mixture  fractions  away  from  the  mean. 

The  noted  inertial-particle  behaviour  is  even  stronger  in  nature  than  that  described  in 
Section  4.1.1  for  smaller  soot  particles.  This  suggests  that  the  larger  the  particles  under 
consideration,  the  less  effective  a  continuum  model  will  be  in  describing  their  behaviour. 


5  Remarks 


The  results  of  this  study  indicate  that  a  tailored  variant  of  the  Conditional  Moment 
Closure  (CMC)  method  can  be  applied  successfully  to  modelling  the  evolution  of  soot¬ 
like  and  droplet-like  reacting  particles  in  a  turbulent  combusting  environment.  For  the 
soot  particle  cases,  a  single  set  of  constants  were  sufficient  to  correctly  predict  particle 
evolution  for  a  variety  of  reaction  intensities.  This  finding  suggests  that  it  is  appropriate 
to  simultaneously  model  condensed  and  gas  phase  reactions  in  turbulence  using  the  CMC 
model  within  an  Eulerian  frame. 

However,  significant  caveats  apply  to  this  modelling  approach.  Firstly,  it  is  clear 
that  the  conditional  covariance  between  scalar  dissipation  rate  and  reactive  scalar  mass 
fraction  cannot,  in  general,  be  neglected.  It  is  apparent  that  some  species  can  exhibit 
strong  non-zero  covariance  profiles  which  play  the  foremost  role,  of  all  transport  terms, 
in  transporting  scalar  properties  in  mixture  fraction  space.  An  examination  of  simulation 
cases  with  different  particle  reactivity  suggest  that  those  species  subject  to  more  intense 
reactions  exhibit  stronger  covariance  profiles.  This  may  be  due  to  the  establishment  of 
large  changes  in  scalar  values  and  gradients  between  adjacent  mixture  fractions,  which  can 
result  from  strong  reactions  taking  place.  It  was  found  that  changes  in  particle  inertia, 
while  keeping  all  else  constant,  had  virtually  no  effect  on  the  form  of  the  conditional 
covariance  profile.  Further  work  is  obviously  required  to  study  the  dependence  of  mixing 
covariance  on  particle  reactivity,  and  to  devise  an  improved  model  for  predicting  its  profile 
shape  in  mixture  fraction  space. 

Secondly,  it  was  found  that  the  CMC  model  is  not  well  suited  to  predicting  the  be¬ 
haviour  of  particles  with  substantial  levels  of  inertia.  Relatively  large  particles,  such  as  fuel 
droplets,  can  develop  significant  slip  velocities  relative  to  their  surrounding  fluid.  These 
slip  velocities  appear  as  unclosed  transport  and  variance  production  terms  in  the  CMC 
equations.  They  appear  to  enhance  transport  of  particle  properties  in  mixture  fraction 
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Figure  9:  Observed  conditional  mean  and  variance  of  particle  mass  fraction  for  inertial 
(dl)  and  inertia-less  (dli )  droplet  cases.  Symbols  denote  the  following  profiles:  A  -  Qv 
inertial,  o  -  Qv  inertia-less,  o  -  qn  inertial,  and  x  -  qv  inertia-less. 
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space  over  that  exhibited  by  smaller  particles  with  smaller  slip  velocities,  but  similar  re¬ 
activity.  Further,  as  one  would  expect,  because  the  larger  particles  are  not  transferred  in 
space  in  the  same  way  as  gaseous  mixture  fraction,  the  conditional  variance  of  particle 
properties  at  any  given  mixture  fraction  is  higher  than  that  for  gaseous  species.  While  it 
was  found  that  the  modelling  of  inertial  particle  transport  in  mixture  fraction  space  could 
largely  be  accounted  for  by  increasing  the  level  of  diffusive  transport  (through  increasing 
c  to  unrealistically  high  levels),  the  model  could  not  account  for  the  increased  conditional 
variance  which  was  observed. 

The  mixing  mode  variable  (c)  approach  to  treating  differential  diffusion  in  the  CMC 
model  deserves  further  attention.  In  varying  the  value  of  c  between  the  extreme  value 
given  by  the  ratio  of  molecular  diffusivities  and  unity,  according  to  changes  in  global 
mixing  conditions,  the  CMC  equations  change  in  a  natural  way  to  embody  differential 
diffusion  effects  to  a  greater  or  lesser  degree  as  required.  A  single  value  of  c  was  found 
to  suffice  for  the  prediction  of  particle  behaviour  in  the  cases  studied  here  under  uniform 
mixing  conditions,  but  with  widely  different  reactive  behaviour.  This  suggests  c  may  well 
be  independent  of  local  variations  in  mixture  fraction  and  particle  mass  fraction,  and  a 
function  only  of  global  mixing  parameters.  Further  work  is  required  to  properly  test  this 
hypothesis,  and  the  model  for  c  proposed  in  Eqn  43. 

Future  work  involving  the  CMC-particle  model  derived  here  will  focus  on  validation 
of  the  model  against  larger  simulations  which  incorporate  more  realistic  chemistry,  higher 
turbulence  levels,  and  a  more  diverse  set  of  mixing  and  reaction  conditions. 
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